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1 Comparing sex and hermaphroditism 

1 compared sexual(SX) and hermaphroditic(HA) reproduction schemes within 
the Penna model [1]. A fixed set of parameters was used, to make the results 
comparable with recent simulations [2]. I also reproduced the results for a self- 
emerging dominance [2] in an independant simulation for both hermaphroditic 
and sexual cases. A variation on partner selection was also implemented, this 
shifted stable population sizes slightly. 

The stable population is considered the most significant result for these com- 
parisons since Stauffer et al.[3] gave strong argument, that from two species com- 
peting for the same resources the one with a higher stable population while not 
competing will prevail, while the other one will die out. 

The parameters were: Birthrate B = 4, minimum reproduction age R = 8, 
mutational threshold T = 3, and mutation rate M — 1. Different values for N max 
have been simulated. The starting population was choosen to N(0) = N max /100. 

In order to get better statistics several runs have been made with a different 
seed value for the random number generator, when computer time was available. 

2 Constant dominance 

To take the effect of dominant and recessive genes into account, a dominance 32 
bit long string was generated for the whole population with d = 6 randomly set 
bits. Deletrious mutations only take effect if either both genome strings have the 
according bit set or one string has it set and the bit is marked as dominant by 
the dominance string. 
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Figure 1: Comparing sexual (top curve) and hermaphroditic (bottom curve) 
reproduction schemes. The dominace was assumed to be constant for the whole 
population for these simulations. 



As can be seen in table EJthe hermaphroditic populations consistantly reached 
a stable population of N Q /N max w 0.23 whereas the sexually reproducing popu- 
lations reached N /N max fh 0.17. No significant difference was found for longer 
simulation time or larger populations. 
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Figure 2: Comparing sexual(SX) and hermaphroditic (HA) reproduction 
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Figure 3: (x) is the number individuals with a dominance bit set at that age, 
(+) is the total number of individuals at that age (H, J, N ~ 0.23) 



3 Dominance self emergence 

I started from the standard Penna model with sexual and hermaphroditic repro- 
duction. Each individual now has a dominance bit string attached. With every 
reproduction one bit in this string is flipped with a probability of p — 0.01. This 
dominance string is then used to determine the activity of a deletrious mutation 
just as in the constant dominance model. As in recent simulations [2] the number 
dominant bits nd om reached a stable value after t > 10 5 time steps (figure 6) and 
has an age distribution reflecting the selection pressure gradient (figure El). 

3.1 Partner selection style 

In recent simulations [3] each female or hermaphroditic individual with age > R 
trys to find a fitting partner from the whole population 20 times. If no partner is 
found, the individual chooses to not reproduce that time step. The simulations 
with this reproduction style are marked 20T. 
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Figure 4: From top to bottom: a.) Hermaphroditic reproduction, J c.) 
Hermaphroditic reproduction, 20T b.) Sexual reproduction, J d.) Sexual re- 
production, 20T 
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Figure 5: Dominance emergance simulation parameters and results 
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Figure 6: Number of dominant bits averaged over the whole population. The 
fluctuating curve represents one run, whereas the smooth curve results from av- 
eraging over 50 runs. (SX,20T) 

Another model was simulated, where each turn a pool of males in the apropi- 
ate reproduction age is generated. Each female chooses one from this pool and 
removes it from this pool. So the individuals form couples each turn. This model 
is marked as J according to a german saying that each jacket will find fitting 
trousers. 

However this different partner selection style only increased the stable popu- 
lation sizes slightly. 

3.2 Population 

As can be seen in figure 4 in the hermaphroditic case a stable population of 
No/N max ps 0.23 was reached, whereas the sexual reproduction scheme reached 
N /N max ~ 0.17, The different partner selection styles simulated did not produce 
important differences. 
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4 Summary 



All simulations consistantly showed that hermaphroditism reaches a higher stable 
population size than sexual reproduction. The Penna model thus does not explain 
the evidence found in nature that sexual reproduction offers severe benefits over 
hermaphroditic reproduction. This known feature of the Penna model has been 
confirmed. 

Several researchers have tried to give an explanation for sexual reproduction in 
comparison to asexual reproduction. Other reproduction schemes (especially mei- 
otic parthenogenesis) have also been analysed. See chapter 2.5.5 in Evolution, 
Money, War and Computers by Oliviera, Oliviera and Stauffer [5]. 

To further compare reproduction schemes one should investigate the applica- 
bility of these approaches to hermaphroditism. 
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